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ANALYTICAL  CHARACTERIZATION  OF  BISTATIC 
SCATTERING  FROM  ROUGH  SURFACES:  DEPENDENCE 

ON  SURFACE  CORRELATION  FUNCTION 


1.  INTRODUCTION 


The  short  pulse  width  In  high  resolution  radars  has  the  effect  of  shortening  the  clutter  cell 
Illuminated  by  the  radar.  This  smaller  area  would  result  tn  smaller  clutter  returns  If  the 
statistical  distribution  of  the  power  scattered  by  the  clutter  does  not  change  as  the  cell  size  is 
decreased.  The  scattered  power  from  large  clutter  cells,  which  have  many  similar  sized 
scattering  facets,  has  a  Rayleigh  distribution.  A  large  number  of  scattering  facets  Is  a 
requirement  for  Rayleigh  scattering.  However,  as  the  cell  size  is  decreased,  the  number  of 
scattering  facets  Is  decreased,  and  Rayleigh  scattering  may  no  longer  result.  In  this  report,  we 
will  see  what  effect  decreasing  the  clutter  cell  size  has  on  several  different  types  of  surfaces 
In  studies  by  Papa  and  Woodworth1  and  Sharpe2  it  is  shown  that  decreasing  the  size  of  the 
clutter  cell  can  significantly  affect  the  statistical  distribution  of  the  scattered  power.  These 
studies  have  only  examined  surfaces  which  are  described  by  a  Gaussian  correlation  function. 
In  this  report  other  surfaces  that  are  described  by  a  power  law  correlation  function  and  an 

Received  for  Publication  12  June  1992 

1  Papa,  R.J.  and  Woodworth,  M.B.  (1991)  The  Mean  and  Variance  of  Diffuse  Scattered  Power 
as  a  Function  of  Clutter  Resolution  Cell  Size.  RADC-TR-91-09. 

2  Sharpe,  L.M.  (1991)  Analytical  Characterization  of  Bistatic  Scattering  from  Gaussian 
Distributed  Surfaces.  RL-TR-9 1-351. 
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exponential  correlation  function  are  studied.  It  Is  found  that  the  type  of  surface,  whether  It  Is 
smoothly  varying  or  very  jagged,  has  a  significant  impact  on  the  distribution  of  the  scattered 
power. 

The  first  section  of  this  report  will  give  a  brief  description  of  the  correlation  functions  that 
are  used,  as  well  as  the  type  of  surfaces  described  by  each  one.  Results  for  the  different 
correlation  functions  will  then  be  presented  and  compared.  Finally,  some  conclusions  about 
the  effect  of  cell  size  and  correlation  functions  on  clutter  statistics  will  be  made. 


2.  CORRELATION  FUNCTIONS 


A  complete  statistical  description  of  a  random  rough  surface  Involves  the  correlation 
function  of  the  surface  as  well  as  the  height  distributions.  The  height  distribution  functions 
at  each  point,  which  describe  the  height  deviations  from  a  certain  mean  level,  are  assumed 
here  to  be  Gaussian  with  a  mean  level  of  zero.  The  correlation  function,  C(t),  is  a  function  of 
the  separation  of  two  points  on  the  surface  and  is  a  measure  of  how  closely  related  the 
distributions  at  these  two  points  are.  The  separation  parameter,  t.  Is  given  by  t  =  (xj  - x 2  ).  with 
and  X2  being  points  on  the  surface.  When  C(t)s  1.  the  heights,  ^  and  C,2  at  xi  and  x2  are 
fairly  strongly  correlated  so  that  knowledge  of  the  height  at  one  point  gives  information  about 
the  height  at  the  other  point,  and  when  C(t)s  0.  the  heights  are  uncorrelated,  so  that 
knowledge  of  the  height  at  one  point  tells  nothing  of  the  height  at  the  other  point.  For  a 
purely  random  surface.  C(t)  is  an  even  function  and: 


Mnl  C(T)  =  1.  l,m  C(t)  =0  . 
T->0  T-><» 


The  correlation  distance.  T.  Is  defined  to  be  the  distance  In  which  C(t)  drops  to  the  value  e-1.  It 
is  this  correlation  distance  that  determines  whether  a  clutter  cell  is  large  or  small.  For 
inllnite  clutter  cells.  T  <<  L  where  L  Is  the  dimension  over  which  the  scattering  occurs.  The 
total  scattered  power  divided  by  L  Is  the  average  scattered  power.  However,  tn  the  finite  cell 
size  studies  done  here,  the  cell  size  ranges  lrom  L=2T  to  L=12T,  so  the  conditions  necessary  for 
Rayleigh  scattering  are  violated  to  different  extents. 

The  correlation  functions  used  in  this  study  are  the  Gaussian.  Power  Law.  and  Exponential. 
They  are  given  by: 


Cc(t)  =  exp  M2/T2)  Gaussian 


(2) 


2 


Cp,(T)  =  [1  +  (t/T)2)-3/2 


Power  Law 


(3) 


Cexp(r)  =  exp  (—1x1  /T)  Exponential 


14) 


The  Gaussian  correlation  is  used  to  describe  smoothly  curving  surfaces  having  height 
derivatives  at  all  points.3  The  power  law  correlation  represents  surfaces  that  are  somewhat 
less  correlated  than  the  Gaussian  surface.  Asphalt  roads  and  gravel  surfaces  are  typically 
described  by  a  power  law  correlation  function.4  Barrick  and  Peake5  found  a  power  law 
correlation  empirically  fitted  measured  points  of  an  asphalt  surface.  Surfaces  described  by  an 
exponential  correlation  function  are  jagged  and  have  many  vertical  facets,  such  as  urban  areas 
having  buildings  and  houses  as  the  scattering  surfaces.  The  three  different  correlation 
functions  are  plotted  in  Figure  1.  The  correlation  distance  is  T=1  for  all  three  cases.  Note  that 
at  t=l,  all  three  functions  cross  the  e-1  point.  Surfaces  described  by  these  different  correlation 
functions  are  quite  different  and  can  have  quite  different  scattering  statistics. 

Figures  2  through  7  show  two  different  realizations  for  surfaces  with  the  three  correlation 
functions  studied.  The  rms  height  of  these  surfaces  is  o=0.5m,  which  is  obviously  not  a 
realistic  value  for  asphalt  surfaces,  where  one  would  expect  a  «  0.5m.  or  urban  areas  where 
o  »  0.5m.  but  for  comparing  the  different  correlation  functions,  a  common  value  for  n  was 
used.  To  determine  these  sample  surfaces  a  random  set  ol  uncorrelated  Gaussian  heights  was 
generated  and  multiplied  by  the  appropriate  correlation  matrix.  The  generation  of  the 
correlation  matrix  and  correlated  surfaces  is  discussed  in  Appendix  A.  The  result  will  be  the 
correlated  surface  heights.  For  Figures  2,  3  and  4.  the  same  set  of  random  numbers  was  used  to 
generate  the  surfaces.  In  comparing  these  three  figures  it  is  obvious  that  the  exponential 
surface  is  much  more  Jagged  and  has  many  more  scattering  facets  than  the  Gaussian  correlated 
surface.  The  surface  with  a  power  law  correlation  is  not  as  Jagged  as  the  exponential  and  not 
as  smooth  as  the  Gaussian  surface.  Even  though  the  surfaces  are  all  quite  different,  the  same 
large  scale  structure  can  be  seen  in  all  three  because  the  same  set  of  random  numbers  was  used 
to  generate  all  tliree  surfaces.  The  surfaces  shown  in  Figures  5,  6  and  7  were  generated  by 
different  sets  of  random  numbers.  Again  we  see  that  the  Gaussian  correlated  surface  is  fairly 
smooth,  the  exponential  surface  Is  quite  Jagged,  and  the  power  law  correlated  surface  lies 
somewhere  between  the  other  two  in  terms  of  irregularities.  Because  these  three  surfaces  were 
generated  by  different  sets  of  random  numbers,  the  large  scale  structure  is  not  common  to  the 
three  surfaces. 


3  Ruck.  G.T.,  Barrick,  D.E.,  Stuart.  W.D.,  and  Krichbaum.  C.K.  (1970)  Radar  Cross  Section 
Handbook.  2,  Plenum  Press.  New  York. 

4  Cohen.  E.,  Bell,  M.,  Un,  C.,  Holst.  J.,  and  Tomasanis.  D.  (1988)  Mathematical  Analysis  and 
Programming  for  Electromagnetic  Applications.  RADC  -TR-  88-18.  ADB 1 24496 . 

5  Barrick,  D.E.,  and  Peake,  W.H.  (1967)  Scattering  from  Surfaces  with  Different  Roughness 
Scales:  Analysis  and  Interpretation,  BAT- 197A- 10-3. 
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Figure  2.  Integration  Area  Being  Filled  by  a  Sobol  Sequence 
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Figure  4.  Power  Law  Correlation  Surface 
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Figure  7.  Power  Law  Correlated  Surface 


3.  ONE-DIMENSIONAL  SCATTERING  MODEL 


Scattering  from  large  clutter  cells  is  typically  assumed  to  have  a  Rayleigh  distribution. 
The  equation  for  a  Rayleigh  distribution  is: 

p<x)  =  1  exp  (-x/g) , 


(X)2,  where  x  is 


where  x  is  the  power.  The  mean  is  g.  and  the  variance  is  given  as  a00  =  (x2 j  - 
the  mean,  and  (x2^is  the  second  moment  given  by: 


oo 


X  exp  (-  x  dx  =  2g2  . 
M  V  gj 


(5) 


The  variance  is  therefore  given  as  ct°°=  2g2  -  g2  =  g2.  and  the  relationship  between  the  mean 
and  the  variance  of  the  scattered  power  for  Rayleigh  scattering  is  o00  =  (o0)2.  However,  as  die 
cell  size  is  decreased,  violating  the  condition  T«  L,  this  relationship  may  no  longer  be  valid. 
In  previous  studies,  the  mean  value  of  the  scattered  power  for  a  fintte  clutter  cell  of  a  one 
dimensionally  rough  surface  was  given  as1-2: 


o 

o  = 


n  F2 
2LA. 


!,/2  L/2 

j  ClXl  j  <^2 

-L/2  -L/2 


*2 -*1*1  cos(vxt) 


(6) 


where: 


X  =  Xj  -  X2 

X2  =  exp|-L2  (1-C(t)|  is  (he  characteristic  function  for  bivariate  height  distribution 

X 1  =  exp  (-£2/2)  is  the  characteristic  function  for  univariate  height  distribution 

2 

X2  =  o2vz  =  Rayleigh  parameter  squared 

0  =  rms  surface  height 

vx  =  (2/r/X)[sin(0,)-sin(03)| 

vz  =  -(2x/X.)lcosl0,)+cos(09)l 


2R  (1  +  cos  +  4>g)) 
cos  0,  +  cos  6 

1  s 


However,  performing  a  change  of  variables  from  x,  and  to  t  and  o  makes  it  possible  to 
solve  (he  integral  over  a  analytically,  leaving  a  single  Integral  for  the  mean  value  of  the 
scattered  power  given  by  (See  Appendix  B): 


o 

o  = 


K  F 


LX 


L 

Jdt[b-i  X2  -  Z,xl  «>s(vxt). 


0 


(7) 


The  variance  for  the  finite  clutter  cell  is  calculated  from: 


L/2 

lllldXlCtX2dX,dX4  008  (Vx  (X1-X2+X3- X4»  +  f*4  "  *12  *34 

\n 

where 


00 
o  =  i 


n  F 
2LX 


(8) 


X12  =  exp  [-I2  ( 1  — C ,  2)1 

X34  =  exp  (-X2  (i-c:  M)l 

Ct)  =  C(t)  as  defined  in  Eq.  (2).  (3),  or  (4). 


The  calculated  variance  is  compared  to  the  Rayleigh  value  of  the  variance,  000  =  (00)2  to 
determine  if  decreasing  the  cell  size  will  cause  the  scattering  to  become  non-Rayleigh.  Note 
•hat  the  surface  correlation  function  is  used  in  the  calculation  of  both  the  mean  and  the 
variance  ot  the  scattered  power. 
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Because  of  the  Intensive  calculations  done  in  this  study,  particularly  for  the  variance, 
efficient  numerical  integration  techniques  were  required.  The  Integration  scheme  used  is  a 
Sobol  sequence.6  It  is  considerably  more  efficient  than  conventional  integration  schemes.  In 
this  method,  the  n-space  that  is  being  integrated  over  is  filled  by  quasi-random  points  that  are 
"maximally  avoiding”  each  other,  and  therefore  fill  the  space  more  uniformly  than  a  purely 
random  sequence.  Figure  8  shows  a  two-dimensional  space  being  filled  by  Sobol  sequences 
The  algorithm  used  for  Sobol  sequence  integration  as  well  as  a  description  of  the  technique  is 
given  in  Reference  5. 


4.  GAUSSIAN  AND  POWER  LAW  CORRELATION  FUNCTION  RESULTS 


The  normalized  mean  and  variance  of  the  scattered  power  for  the  Gaussian  correlated 
surfaces  are  shown  In  Figures  9  through  16.  The  curve  with  the  small  boxes  represents  the 
mean  value  of  the  scattering  a0,  the  curve  with  the  triangles  represents  the  Rayleigh  value  ot 
the  variance  (o°)2,  and  the  curve  with  the  large  boxes  represents  the  calculated  variance  for  the 
finite  clutter  cell.  In  Figures  9.  10.  and  1 1.  0,=75°,  A.=0.25m.  o=0.5.  and  L  decreases  from  12m 
in  Figure  9,  to  fim  in  Figure  10.  and  to  3m  in  Figure  11.  It  is  clear  that  as  the  cell  size 
decreases,  the  scattering  becomes  less  Rayleigh,  particularly  in  the  backscatter  region,  due  to 
fewer  scattering  facets  oriented  in  this  direction.  This  means  that  decreasing  the  cell  size  for 
this  type  of  surface  will  have  more  severe  consequences  for  a  monostatic  radar  than  for  a 
bistatic  radar. 

Figure  12  is  the  same  as  Figure  11  except  that  the  wavelength  is  decreased  to  >.=0.1  m 
Comparing  Figures  1 1  and  12  shows  that  the  scattering  is  closer  to  Rayleigh  for  the  longer 
wavelength.  This  is  because  there  are  more  scatterers  per  wavelength  as  the  wavelength  is 
increased. 

Figure  13  is  the  same  as  Figure  12  except  the  surface  roughness  is  decreased  to  o°-0.2m. 
For  the  smoother  surface  the  scattering  is  concentrated  more  in  the  specular  direction  while 
the  rougher  surface  scatters  more  diffusely. 

For  Figures  14  through  16  the  angle  of  incidence  is  01=30°.  These  results  show  the  same 
trends  as  seen  for  the  8,=75°  case,  that  is,  for  the  longer  wavelength  the  scattering  is  (  loser  to 
Rayleigh,  and  the  scattering  from  the  rougher  surface  is  more  diffuse  over  the  scattering 
region. 


6  Press,  W.H.  and  Teukolsky.  S.A.  (1989)  Quasi-  (that  Is  sub  -)  random  numbers.  Computers  in 
Physics:  76-79. 
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Figure  8.  Exponential  Correlated  Surface 
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o  =  0.5m,  L  =  12m 
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SCATTERING  ANGLE  (degrees) 
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Figure  1 1.  Normalized  Mean  and  Variance  for  Gaussian  Surface.  0,  =  75°.  X  =  0.25m, 
o  =  0.5m.  L  =  3in 
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SCATTERING  ANGLE  (degrees) 

Figure  13.  Normalized  Mean  and  Variance  for  Gaussian  Surface.  0,  =  75°. 
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Figure  14.  Normalized  Mean  and  Variance  for  Gaussian  Surface,  9j  =  30°,  k  =  0.25m, 
a  =  0.5m.  L  =  3m 
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SCATTERING  ANGLE  (degrees) 

Figure  15.  Normalized  Mean  and  Variance  for  Gaussian  Surface,  0t  =  30c 
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Figure  16.  Normalized  Mean  and  Variance  for  Gaussian  Surface,  0t  =  30°,  X  =  0.25m, 
o  =  0.2m,  L  =  3m 


The  scattering  from  surfaces  having  a  power  law  correlation  function  are  shown  in 
Figures  17  through  24.  The  results  for  the  power  law  correlated  surfaces  are  quite  similar  to 
those  obtained  with  a  Gaussian  correlated  surface.  The  trends  due  to  changing  wavelength, 
rms  height,  cell  size,  and  scattering  angles  are  all  the  same  for  surfaces  with  a  power  law 
correlation  function  as  they  were  for  the  Gaussian  correlated  surface.  One  notable  difference 
can  be  seen  In  comparing  the  results  for  the  two  surfaces.  It  is  clear  that  in  the  backscatter 
region,  where  scattering  from  both  surfaces  deviates  from  Rayleigh,  the  deviation  for  the 
power  law  case  is  not  as  severe.  This  Is  a  consequence  of  the  Increased  number  of  scattering 
facets  oriented  to  cause  backscatter  for  the  surface  with  a  power  law  correlation.  In  comparing 
Figures  3  and  6  to  Figures  2  and  5,  it  is  obvious  that  the  power  law  correlated  surface  is 
somewhat  more  Jagged  and  will  therefore  have  more  scattering  facets  than  the  Gaussian 
correlated  surfaces.  The  increase  In  the  number  of  scattering  facets  results  in  scattering 
statistics  that  are  closer  to  Rayleigh. 


5.  EXPONENTIAL  CORRELATION  FUNCTION  RESULTS 


The  surfaces  described  by  an  exponential  correlation  function,  as  stated  previously,  are 
very  Jagged  and  have  many  vertical  facets.  Also,  this  type  of  surface  has  more  fine  scale 
structure  than  the  other  surfaces.  Because  of  this,  the  scattering  from  this  type  of  surface  is 
quite  different  than  that  of  the  Gaussian  or  power  law  correlated  surfaces.  The  most 
significant  result  found  for  the  exponentially  correlated  surface  Is  that  the  cell  size  has  very 
little  effect  on  the  Raylelghness  of  the  scattering  statistics.  The  results  for  the  exponentially 
correlated  surface  are  shown  tn  Figures  25  through  32.  For  cell  sizes  of  L  =  3m.  no  noticeable 
deviation  from  Rayleigh  scattering  occurs  over  the  entire  scattering  region.  Because  at  L  =  3m 
the  agreement  with  Rayleigh  scattering  was  so  close,  larger  cell  sizes  were  not  examined.  The 
smallest  cell  size  studied  Is  L  =  2m.  In  this  case  we  see  that  there  is  still  very  little  deviation 
from  Rayleigh  scattering.  This  is  a  very  significant  result  because  it  means  that  for  very  small 
clutter  cells.  Rayleigh  scattering  can  be  assumed  even  for  monostatic  radars  when  the  surface 
is  described  by  an  expone’.-.d  correlation  function.  In  the  Gaussian  and  power  law  cases,  the 
deviation  from  Rayleigh  scattering  is  quite  large  in  the  backscatter  region. 
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SCATTERING  ANGLE  (degrees) 

Figure  18.  Normalized  Mean  and  Variance  for  Power  Law  Surface,  0j  =  ' 
o  =  0.5m,  L  =  6m 
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SCATTERING  ANGLE  (degrees) 

Figure  20.  Normalized  Mean  and  Variance  for  Power  Law  Surface.  0, 
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SCATTERING  ANGLE  (degrees) 

Figure  21 .  Normalized  Mean  and  Variance  for  Power  Law  Surface,  0,  = 
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SCATTERING  ANGLE  (degrees) 

Figure  23.  Normalized  Mean  and  Variance  for  Power  Law  Surface.  0,  *  3 
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SCATTERING  ANGLE  (degrees) 

Figure  24.  Normalized  Mean  and  Variance  for  Power  Law  Surface.  0,  =  ; 
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SCATTERING  ANGLE  (degrees) 

Figure  25.  Normalized  Mean  and  Variance  for  Exponential  Surface.  0,  =  75°.  \  =  0.25m, 
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SCATTERING  ANGLE  (degrees) 

Figure  26.  Normalized  Mean  and  Variance  for  Exponential  Surface.  0,  =  1 
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Mean  and  Variance  for  Exponential  Surface,  9j  =  7 
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Figure  28.  Normalized  Mean  and  Variance  for  Exponential  Surface,  0j  = 


NORMALIZED  MEAN  AND  VARIANCE  OF  SCATTERED  POWER 

EXPONENTIAL  CORRELATION 


IT) 


(ap)  aoNviavA  qnv  nvbiai 


I 

N 

6 

ii 


O 

CO 


36 


SCATTERING  ANGLE  (degrees) 

Figure  29.  Normalized  Mean  and  Variance  for  Exponential  Surface,  0,  = 
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SCATTERING  ANGLE  (degrees) 

Figure  30.  Normalized  Mean  and  Variance  for  Exponential  Surface.  Bj 


SCATTERING  ANGLE  (degrees) 

Figure  31.  Normalized  Mean  and  Variance  for  Exponential  Surface,  8,  =  30°,  X  =  0.25m, 
a  =  0.2m,  L  =  3m 
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SCATTERING  ANGLE  (degrees) 

Figure  32.  Normalized  Mean  and  Variance  for  Exponential  Surface,  0,  =  30°,  k  =  0.25m. 
o  =  0.2m,  L  =  2m 


The  reason  why  the  exponentially  correlated  surface  exhibits  Rayleigh  scattering  for  small 
cell  sizes  while  the  other  surfaces  do  not  can  be  understood  by  considering  the  geometry  of  the 
exponentially  correlated  surface  shown  in  Figures  4  and  7.  In  the  Gaussian  and  power  law 
cases,  where  the  surface  is  smoothly  varying,  the  deviation  from  Rayleigh  could  be  quite 
significant  for  small  cell  sizes,  particularly  in  the  backscatter  region,  because  there  are  so  few 
scattering  facets  oriented  properly  to  cause  backscatter.  However,  the  exponential  surface, 
which  is  very  Jagged  and  has  many  vertical  facets,  has  a  sufficient  number  of  scattering 
surfaces  oriented  in  all  directions  to  cause  the  scattering  to  remain  Rayleigh  over  the  entire 
region  of  scattering  angles.  For  the  exponential  surface,  the  direction  In  which  the  maximum 
scattering  occurs  Is  -90°,  due  to  the  large  number  of  vertical  facets,  while  the  minimum 
scattering  is  in  the  specular  direction. 

The  mean  value  of  the  scattering  from  an  exponentially  correlated  surface  is  almost 
always  lower  than  o°  for  the  Gaussian  correlated  surface.  For  the  30°  angle  of  incidence,  when 
the  wave  is  incident  upon  vertical  facets,  a  considerable  amount  of  the  energy  will  be  reflected 
downward,  and  since  second  reflections  are  not  accounted  for  in  this  study,  this  energy 
disappears.  However,  for  the  75°  angle  of  incidence,  much  of  the  scattering  is  reflected  back  by 
the  vertical  facets  causing  very  high  backscatter  and  very  little  forward  scatter.  In  comparing 
the  mean  value  of  the  scattering  from  the  Gaussian  correlated  surface  and  the  exponential 
surface  we  see  that  o°  is  higher  (or  the  exponential  correlation  In  the  backscatter  region  but 
lower  than  scattering  from  Gaussian  correlated  surfaces  in  the  forward  scatter  region.  This 
result  agrees  with  intuition  because  higher  backscatter  would  be  expected  from  urban  areas  for 
low  grazing  angles.  Table  1  compares  o°  for  the  backscatter  direction  for  the  exponential  and 
Gaussian  correlated  surfaces. 


Table  1.  Comparison  of  Backscatter  for  Exponential 
and  Gaussian  Correlated  Surfaces. 


of 

A.  (m) 

o  (m) 

L  (m) 

°L  <dB) 

exp 

a  (dB> 

ftauss 

30 

0.25 

0.5 

3 

-8.50 

2.90 

30 

0.1 

0.5 

3 

-12.48 

2.95 

75 

0.25 

0.5 

3 

8.79 

-9.18 

75 

0.1 

0.2 

3 

6.81 

-13.65 

6.  CONCLUSIONS 


In  this  report  It  is  shown  that  as  the  size  of  a  clutter  cell  is  decreased,  the  distribution  of 
the  power  scattered  from  that  clutter  cell  depends  on  several  factors.  The  wavelength  of  the 
incident  field,  the  roughness  of  the  surface,  and  the  incident  angle  will  all  have  some  effect  on 
the  Ravlelghness  of  the  scattered  power,  as  discussed  In  the  report  by  Sharpe2,  but  the  results 
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shown  in  this  report  demonstrate  that  the  correlation  function  of  the  surface  can  have  a 
dramatic  Influence  on  the  statistics  of  the  scattered  power.  For  surfaces  described  by  a 
Gaussian  correlation  function,  deviation  from  Rayleigh  scattering  can  be  significant, 
particularly  in  the  backscatter  region.  For  surfaces  having  a  power  law  correlation  function, 
the  scattering  also  shows  a  noticeable  departure  from  Rayleigh  in  the  backscatter  region, 
although  the  difference  is  not  quite  as  large  as  for  the  Gaussian  correlated  surfaces 

For  surfaces  having  an  exponential  correlation  function,  the  scattering  statistics  did  not 
deviate  from  Rayleigh  over  the  entire  scattering  region  studied,  even  for  very  small  sizes.  This 
is  not  surprising  because  there  is  still  a  large  number  of  scattering  facets  in  the  exponentially 
distributed  surface,  even  for  the  very  small  cell  sizes,  which  causes  the  scattering  to  remain 
Rayleigh,  while  the  number  of  scattering  facets  on  the  Gaussian  distributed  surface  decreases 
rapidly  with  cell  size.  From  these  results  we  can  conclude  that  the  number  of  scattering  facets 
in  a  clutter  cell  has  a  significant  effect  on  the  statistics  of  the  scattered  power. 

Knowledge  of  the  type  of  terrain  which  is  to  be  encountered  in  a  radar  system  can  therefore 
be  very  important.  If  the  expected  clutter  is  known  to  have  a  Gaussian  correlation  or  power 
law  correlation,  problems  may  be  encountered  due  to  large  returns  from  the  clutter.  However, 
if  the  surface  causing  the  clutter  Is  known  to  have  an  exponential  correlation  function,  the 
assumption  of  Rayleigh  scattering  will  be  valid,  even  for  very  small  cell  sizes.  Because  an 
exponential  correlation  function  is  used  to  describe  cities  and  other  urban  areas,  this  is  an 
important  result. 
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Generation  of 


the  Correlation 


Matrix  and 


Appendix  A 

Correlated  Surfaces 


To  generate  a  surface  having  a  Gaussian  height  distribution  function  and  a  specified 
correlation  function,  C(x).  one  must  take  the  desired  correlation  matrix,  transform  it  into  an 
upper  triangular  matrix,  and  multiply  it  by  the  set  of  Independent  Gaussian  distributed 
heights.*1  Let  u  =  |u,,  u2,  .  .  .  u^  be  the  set  of  uncorrelated,  Gaussian  distributed  random 
variables,  and  C  be  the  correlation  matrix. 

The  correlation  matrix  is  derived  directly  from  the  correlation  function.  C(t).  The 
demlnsions  of  the  matrix  are  equal  to  the  number  of  points,  N,  on  the  surface  of  length.  L.  The 
distance  between  two  adjacent  points  on  the  surface  is  d=L/(N-l).  t  is  equal  to  the  distance 
between  point  i  and  point  J  and  is  given  by  x=d  I  t-j  I .  The  matrix  can  then  be  filled  in  by  using 
the  correlation  function  and  making  Cjj=C(x).  The  matrix  will  be  square  and  symmetric. 

For  example,  if  a  surface  is  10  units  in  length  and  we  want  to  determine  the  correlation 
matrix  of  six  points  on  the  surface,  we  would  have  L=10,  and  N=6.  Therefore.  d=2  and  x=2 1  i-J  I . 
For  an  exponential  correlation  function,  C(x)=exp(- 1  x  I  /T)  where  T  is  the  correlation  length 
which  we  will  set  to  one  for  this  example.  We  can  then  fill  in  the  correlation  matrix  by  the 
equation.  C,j=exp(-2 1  i-J  I ).  For  i=J,  C,j=l.  Also,  C,j=Cj,. 


A1  Geist,  J.M.  (1979)  Computer  generation  of  correlated  Gaussian  random  variables.  Proc. 
IEEE ,  67  (No.  5). 
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I i-j I  =  Ij-il  =  1,  CtJ  =  exp(-2)  =  0.135 
1 1-J I  =  Ij-i  I  =  2.  C,j  =  expM)  =  0.018 
l  H I  =  IJ-t  I  =  3.  CfJ  =  exp{-6)  =  0.002 
I  i-j  I  =  Ij-i  |  =  4.  C,j  =  exp(-8)  =  0.0003 
I  i-j  I  =  Ij-i  I  =  5.  Cy  =  exp(-lO)  ~  0 

The  correlation  matrix  for  the  exponential  correlation  function  is  then: 


1 

0.135 

0.018 

0.002 

0.0003 

0 

0.135 

1 

0.135 

0.018 

0.002 

0.0003 

0.018 

0.135 

1 

0.135 

0.018 

0.002 

0.002 

0.018 

0.135 

1 

0.135 

0.018 

0.0003 

0.002 

0.018 

0.135 

1 

0.135 

0 

0.0003 

0.002 

0.018 

0.135 

1 

The  next  step  Is  to  convert  the  matrix  into  upper  triangular  form.  The  method  used  to 
transform  the  correlation  matrix  Into  upper  triangular  form  is  the  Cholesky  decomposition  A2 
Let  A  be  the  upper  triangular  representation  of  the  corelation  matrix.  Subroutine  CHFAC. 
which  is  part  of  the  1MSL  library,  performs  the  decomposition .A3  Subroutine  RNMVN,  also 
from  the  1MSL  library,  generates  psuedorandom  Gaussian  distributed  numbers,  u.  to  simulate 
heights  having  a  Gaussian  distribution  function.  These  heights  are  then  multiplied  by  the 
upper  triangular  matrix  found  in  CHFAC.  resulting  in  a  vector,  y 


y  *  A  x 


(A2) 


Tlie  vector  y  represents  the  heights  of  a  surface  that  is  Gaussian  distributed  and  has  a 
correlation  function.  C(r). 


A2  Stark,  H.  and  Woods.  J.W.  (1979)  Probability,  Random  Processes,  and  Estimation  Theory 
for  Engineers ,  Prentice  Hall.  Inc..  Englewood  Cliffs.  New  Jersey. 

A3  MSL,  Inc.  (1987)  Subroutine  CHFAC  and  Subroutine  RNMVN,  Fortran  Subroutines  for 
Mathematical  Applications,  IMSL,  Inc.,  Houston,  Texas. 
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Appendix  B 

Solving  the  Integral  for  Scattered  Power  Over  a 


The  mean  value  of  the  scattered  power  from  a  one  dimensionally  rough  surface  of  length  L 
Is  given  by: 


o 

a  = 


L/2  L/2 


&.W 


F(t) 


(Bl) 


-L/2  -L/2 


where: 


x=  xrXa 


F(t)  = 


X2~X,X1 


cos  (vxt) 
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Perform  a  change  of  variables: 


T  =  X,  —  Xj  x,  =  (1/2)  (x  +  o) 

0  =  X[+X2  X2  =  (1/2)  (o  -  T) 

9x  j  dx  j  i 

j  _  |  dt  <kr  _  i 

2 

9t  3o 


Determine  new  limits: 


T 

=  X, 

- 

=  L.(.L|: 

max 

1  .max 

2. min 

2  12/ 

T 

=  X, 

-  Xrt 

_ L  _  L _ 

min 

1  .min 

2, max 

2  2 

a 

=  X, 

+  X0 

=l+l=l 

max 

1  .max 

2. max 

2  2 

a 

=  x, 

+  X„ 

=  -L+(_L 

min 

1  .min 

2. min 

2  l  2 

-L 


(B2) 


(B3) 


The  Integral  Is  now  over  Ffr)  as  shown  in  Figure  Bl. 
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region  i 


region  2 


L  t-axis 


-L 


Figure  Bl.  New  Integration  Area 


Region  1: 


Region  2: 


(L  +  x)  F(i)  dx 


(L-  x)  F(x)  dx 


Region  3: 


L  o  L 


-Ut  o 


(L-x)  F(X)  dx 


Region  4: 


o  O  o 

fdxF(x)  f  do  =  f  (L  -  x)  F(X)  dx 


Combining  regions  1  and  4  and  regions  2  and  3  and  including  the  factor  of  (1/2)  from  the 
Jacobian  gives: 


=  |  J  (L  +  t)  F(t)  dx  +  2  J  (L  -  x)  F(t)  dx 

L  o  L 

I  =  J  L  F(t)  dx  +  J  x  F(t)  dt  -  j*  t  F(t)  dt 

-L  -L  ° 


(B4) 


Because  F(x)  is  an  even  function,  we  can  write: 


j*  t  F(t)  dr  =  -J*xF(x)  dt 


and 


(B5) 


L  L 

J*  L  F(X)  dx  =  2J*  LF ( x )  dx  . 

O 


The  total  Integral  can  then  be  written  as: 


L 


F(x )  (L-x)  dx  . 


O 


(B6) 


The  double  integral  for  the  mean  scattered  power  can  then  be  written  as: 


O 

a 


Xa-Xjxjjcos 


(V) 


(L-x)dx  . 


(B7) 
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MISSION 

OF 

ROME  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search i,  development,  test,  and  technology  transition  in  support  of  Air 

3 

Force  Command,  Control,  Communications  and  Intelligence  (C  I )  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 

O 

ESD  elements  to  perform  effective  acquisition  of  C  l  systems.  In  addition, 
Rome  Laboratory's  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DO D  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainability  and  testability. 


